
Polytechnic 

Institute 

(of 


A/a5A -cfi-iBK, on 

AERODYNAMICS 

laboratories 

DEPARTMENT OF MECHANICAL 
and 

AEROSPACE ENGINEERING 


THE APPLICABILITY OF BRILLOUIN SCATTERING 
TO FLOW FIELD DIAGNOSTICS 

NASA-CR-158017 

19790010112 

By J. Laiosa and S. Lederman 
FINAL REPORT 
FEBRUARY 1979 


Grant NSG 1198 ■ 

Approved for public release; 
distribution unlimited. 



POLY M/AE Report No. 79-10 





THE APPLICABILITY OF BRILLOUIN SCATTERING 


TO FLOW FIELD DIAGNOSTICS 


FINAL REPORT 


NASA Technical Library 



by J. Laiosa and S. Lederman 


Prepared For 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Langley Research Center 
Hampton, Virginia 


Grant No. NSG-1198 
S. Lederman, Principal Investigator 


July 1, 1976 - June 30, 1977 ) 


POLYTECHNIC INSTITUTE OF NEW YORK 
Aerodynamics Laboratories 
February 1979 


POLY M/AE Report No. 79-10 





TABLE OF CONTENTS 


Page 


Section 
Nomenclature 
List of Figures 

I - Introduction 

II - Theory 

A. Hydrodynamic limit 

B. Kinetic limit 

III - Experimental Criteria 

IV - Experimental Equipment 

V - Gas Sample Preparation and Technique 

VI - Findings 

VII - Conclusion 

VIII - References 

IX - Bibliography 
Figures 


11 

V 

1 

3 

6 

8 

9 

12 

15 

16 
21 
23 
25 
26 - 


1 



THE APPLICABILITY OF BRILLOUIN SCATTERING 
TO FLOW FIELD DIAGNOSTICS 

FINAL REPORT 
by 

J. Laiosa and S. Lederman 
ABSTRACT 

To fill the void between turbulence theory and experiment; par- 
ticularly in the flow fields consisting of monoatomic gases, 
for example in wind tunnels, means of measuring fluctuating 
quantities are needed. In the area of density fluctuation mea- 
surement, the optical method of Brillouin scattering has been 
suggested. This was based on the theory, that -the Brillouin scat- 
tered intensity is proportional to a function of density. 

This function can be shown to be the density-density correlation 
function, (lyo) p (-k) > . In this investigation the potential of 
this method as a diagnostic tool was studied. Here the density 
fluctuations in gases were sought. Continuous wave lasers and 
interferometers were used as the primary illuminating source and 
scattered light filters respectively. 

NOMENCLATURE 

c speed of light, 2.9973x10^^ cm/sec 

Cp specific heat at constant pressure 

c^ specific heat at constant volume 

Crot' ^vib rotational and vibrational specific heats 

c^^ translational specific heat, for monoatomic gas, c^^=c^ 

D self-diffusion constant 

e dielectric constant 


ij 



/s 


E 







k,k. ,k 

' i' s 
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P 
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r 


R 

S (k,w) 


T,To,Ti 


Fourier-Laplace transform of the dielectric constant 

kinetic distribution function, equilibrium value 

perturbed quantity of the kinetic distribution func- 
tion 

imaginary number, i = 

scattered intensity, incident radiation, central line 
intensity, Brillouin line intensity 

change in wave propagation vector, incident and scat- 
tered wave propagation vector 

Boltzman's constant, 1.38x10 ergs/°K 

molecular mass 

number of molecules 

pressure 

momentum 

kinetic energy 

position vector 

molecular gas constant 

generalized structure factor 

temperature, equilibrium value, fluctuating value, 

T = T +T, 
o 1 
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time 
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rot 
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rotational relaxation time 


vibration relaxation time 

modified rotational relaxation time, r' 


modified vibration relaxation time. 
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speed of sound 

non-dimensional parameter 

general thermodynamic parameter 

polarizability 

thermal diffusivity 

thermal expansion coefficient 

shear viscosity 

bulk viscosity 

classical bulk viscosity 

vibrational bulk viscosity 

rotational bulk viscosity 
density 

wavelength of light 
thermal conductivity 
mean free path 

wavelength of scattered light 
angular frequency of scattered light 

angle between point of observation and electric field 
vector 

shift in angular frequency 
frequency 

ratio of specific heats 

angle between incident wave vector and point of ob- 
servation 
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INTRODUCTION 


In the past probes, such as hot wires, thermo- 
couples, pitot tubes and sampling probes, were used to obtain 
good results for average values of fluid parameters. While it 
is true, some of these methods could be used to obtain the fluc- 
tuating values of fluid parameters, care must be taken to elim- 
inate the fluctuations due to the probe itself. This is not 
easily done since its effects are not always known apriori . 

Therefore one has sought other means of investigating fluid flows. 

This train of thought spurred the development of non-intrusive 
optical techniques^. These methods can be sub-divided into two 
classes, indirect (or "path") methods and direct (or "point") 
methods. In the area a density fluctuation measurement, methods ^ 
of the first type include absorption-emission. Interferometry and 
Holography. The second type include such methods, as Raman 
scattering, Rayleigh scattering and florescence. 

In the "path" methods one has to contend with diagnosing only 
2-dimensional or axis-symmetric flows. This is due to interaction 
of light and fluid occurring along the entire beam length. There- 
fore the v7hole beam is the scattering volume. The techniques for 
extracting information from the scattered light are mathematically 

I 

difficult. Not only are they difficult to implement, these tech- 
niques closely resemble the turbulent modeling. Therefore, why not 

^ 

numerically model the flows instead of experimentally trying to 

measure them. 

In the area of direct or "point" techniques, Raman scattering 
has emerged as a useful and accurate method for determining point- 
wise concentrations in flow fields. But for fluctuation measure- 


1 



ments a costly high power laser is needed. In addition, Raman scat- 
tering is applicable only to polyatomic molecules which are Raman 
active. Monoatomic molecules do not exhibit Raman activity. 

In light of this it would be desirable to have a method that 
would produce a higher level of scattering and at the same time be 
sensitive also to monoatomic molecules. With this, a relatively 
.low cost continuous wave (cw) laser could be used to monitor the 
fluid as a function of time. 

At the outset of this investigation, it was believed Brillouin 

scattering could accomplish this feat. In the past, Brillouin scat- 

2 3 

tering has been used to measure sound speeds in gases ' . 

The procedure used to monitor the density fluctuations involves 
the use of a high resolution interferometer in an unconventio.nal 
mode of operation. Usually a sawtooth is applied to the piezo- 
electric drive which provides the sweep over a narrow bandwidth 
corresponding to the free spectral range of the interferometer. 

If a single voltage is applied, theoretically a single frequency 
should pass through the mirrors. Therefore, by prebiasing the 
interferometer at a particular voltage a particular frequency of 
light will be transmitted. And since the intensity of this light 
is proportional to the density of the fluid, a measure of the 
variation of this scattered light would indicate the fluctuation 
in the density of the fluid. In an 

earlier investigation^ an attempt was made to measure density fluc- 
tuations in liquids, but due to quality of the equipment available, 
and the difficulty of purifying these liquids, no conclusion could 
be drawn at that time. Therefore, new equipment was obtained and 
the interest was shifted to gases since it was believed the gases 
could be easily purified. The interest in gases v/as also sparked 
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by the fact that it is the primary type of fluid used in fluid dynam 
ic flow fields. Liquids were only introduced because their scatter- 
ing is more intense and therefore an easy point to begin an investi- 
gation. 

While a Rayliegh-Brillouin spectrum could not be detected, some 
reasons were discovered that implied the inadequacy of Brillouin 
scattering as a diagnostic tool at the present time. 

In the following section a brief account of the theoretical as- 
pects of Brillouin scattering is shown. The following section pro- 
vides the experimental criteria needed for a successful experiment. 
After a description of the equipment involved, the reasons for dis- 
counting Brillouin scattering to measure density fluctuations are 
given. 
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II - THEORY 


The phenomena of Brillouin scattering is the result of light scat- 

4 

tering from density fluctuations present in the fluid. Brillouin 
first implied that light is elastically scattered from thermal 
sound waves produced by the medium. These sound waves create 
stratified layers consisting of alternating Planes of compression 
and rarefraction wave forms (see Fig. 1) . 

Due to momentum conservation requirements, the change in wave 
propagation vector is given by 


K = K. -K 


1 s 

The magnitude of K can be expressed as 


( 1 ) 


|k1 = 2 |k^ I sin 0/2 


( 2 ) 


since |k^|«|k^1. 

Also, in the presence of the translational motion of the 
molecules, the scattered light is Doppler shifted. This shift xs 
given by 


Av 


^^o 

2tt 


2V sin0/2 
o 

X 


(3) 


o 

The scattered light intensity can be expressed in many ways. 

6 V 

Kamarov and Fisher^ have adapted Van Hove ' s ' neutron-scattering to 
light scattering. Their expression for the scattered light intensity 
is : 


I(r,w) = (a^fi^N/2TTc'^r^) I^san^'I'S 0 ^,( 0 ) 


(4) 
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The generalized structure factor S can be expressed: 

S(k,u) = 2R^<g (k,iu))e (-k) > (5) 

<• 

Where § is the Fourier-space Laplace-time transform of the di- 
electric constant e. Therefore one needs to express the behavior 
of e, for the medium. Fortunately the change in e can be ex- 
pressed as the change in any two thermodynamics parameters. If density 
and temperature are selected, one obtains: 

6e(F,t) = (||)^6p(F,t) + (||) 6T (F,t) (6) 

This can be further simplified since the temperature dependence is 

9 

small as compared to the density dependence . Therefore the prob- 
lem rests on obtaining an expression for 9p(F,t). 

Due to the diverse nature of dense fluids and gaseous fluids, 
two regions exist where different solutions apply. Hence, 3p may 
be obtained in different manners for each region. Before discus- 
sing these methods, it would be advantageous for one. to know in 
which region a particular fluid resides. This is important since a 
gas can be in either region depending on its thermodynamical prop- 
erties and the scattering angle used for viewing. 

To distinguish between the two regions, one should consider the 
collisional rate process which differs vastly for each. With this, 
one can intuitively see the validity of the parameter "y" as an indica- 

g 

tor of the region in which the fluid resides for a particular case . 
This parameter is expressed as: 

y = ^s/^c "" C2iT/k)/X^ (7) 

Here y is seen to be the ratio of wavelength of the fluctuations ob- 
served to the collisional mean free path. 
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Therefore for y>>X the wavplenqth of the. fluctuations is jauch 
larger than the mean free path. This regime is denoted the 
"Hydrodynamic regime" since it is governed by the hydrodynamic 
equations of motion. 

For y<<l/ the fluctuation wavelength is much shorter than 
the mean free path. Since the collisions are very infrequent, the 
system can no longer be considered a continuum. Hence, this is 
denoted as the "kinetic regime". This variation of y can be 
viewed graphically by considering Fig. 2, (Ref. 4). Here the 
parameter y is shown for two monoatomic molecules. Helium and 
Xenon at two scattering angles, 10° and 160°. As stated earlier, 
each molecule can be in either regime depending on k (scattering 
angle) and its density. This can also be seen in the above figure. 

Now that the distinction between the two regimes has been de- 
fined, an expression for 6p for each regime is possible. Before 
these expressions for 6p are given, the spectral shape of the scat- 
tered light for each domain will be discussed. In the hydrody- 
namic region the spectriam consists of three components, (see 
Fig. 3) . The central component is the unshifted Rayleigh line. 
Straddling this central line are the two Brillouin lines. The 
existence of tv;o lines is due to the oscillations of the thermal 
sound v;aves in both longitudinal directions. The ratio of intensi- 
tives of this spectrum can be expressed as 


21 


= Y-1 


B 


( 8 ) 


where I and are the intensities of the central and Brillouin 
c B 

lines respectfully. This equation actually holds where relaxation 
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effects are negligible, but it gives a good result for relative 
differences in the intensities of the central and Brillouin 
peaks. As one proceeds into the kinetic region (y decreasing), 
the lines come closer togetlier. This can be seen in Fig. 4 
taken from Ref. 4 . Here it is seen at a value of y = 2.0 the 
lines are still separated, but at y = 1.0 the lines are almost 
indistinguishable. Finally at y = 0.5, no distinct line is visi- 
ble. Also, the maximum intensity tends to decrease as y decreases. 

In the sections that follow, a brief statement of the deriva- 
tion of 6p for both regimes is given. These derivations are taken 
primarily from Ref. 10 & 11. The interested reader can consult 
these references for a more descriptive explanation. 

A. Hydrodynamic Limit 

For a dense fluid, the linearized hydrodynamic equations of 
12 

motion are the continuity equation: 


3p-| 

^ + Po ° 

the longitudinal component of the Navier-Stokes equation: 

Po gradpjj^ + — VT^ - (-jn^ + n^) grad div u=0 (10) 

and the energy transport equation; 


BT^^ c^(Y-l) 
Po^v 8t p 


8pj^ 

pt“ 


X^V^T^ = 0 


( 11 ) 


From here one can obtain the Four ier-space Laplace-time trans- 
form of the density fluctuations. This can be accomplished by first 
taking the divergence of Eq. 5'' and substituting for V-u^ from 
Eq. 4. Then the Fourier-Laplace transforms can be obtained for 6p 
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(see Refs. 11,13). The final result for the intensity can be 
obtained numerically or by several approximate analytical methods . 

The “above theory provided satisfactory results for monoatomic 
fluids. For polyatomic fluids, the theory differed appreciably 
from experiment. Neglect, of the other forms of internal degrees 
of freedom (rotational and vibrational) was the reason for this 
discrepancy. 

The inclusion of these effects can be made by introducing a 
frequency dependent bulk viscosity and a self diffusion dependent 
thermal conductivity in the given hydrodynamic equations. For the 
bulk viscosity we have 


n 


B 


CLASS 


vib 


+ n 


rot 

B 


( 12 ) 


where n 


B 


CLASS 


is the usual bulk viscosity and 


rot ^^^rot'^^rot 


'B 


(C ^+C^ ) 
rot tr 


2 l+xmx 


(13) 


n 


vib 


PRC T 
vib vib 


B 


1+ltOT 


(14) 


are the rotational and vibrational bulk viscosities respectively' 
As for the thermal conductivity , one has the following expression 


= R [| C. + p D(C . + C ., )] (15) 

T s 2 tr ^o rot vib 

With these modifications, the scattering from a polyatomic fluid 
can be accurately predicted. 

As a final note, theories have been developed for binary gas 
mixtures. These theories can be seen in Refs. 11 and 16. 
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B. Kinetic Limit 


In the kinetic regime, one must express the density distrib- 
ution in a manner compatible to iwith the _theory of rarefied gases. 

For a monoatomic gas, the molecular distribution is given by 

dN(u,r,t) = f(u,r,t)drdu (16) 

Here dN is the number of molecules with the coordinates (u,r,t) 
within the spatial volume dr and velocity voliame du. With this 
distribution, any macroscopic parameter can be determined. This 
can be seen if one considers y(r,u) as any one of those parameters. 
Then 


<y(r, t)> = /f (u,r ,t)y (r,u) du (17) 

If y is a power^o?^ffie^velocity only, it is then considered a moment 
of the function f. In the collisional process there are several 
moments of physical interest. These are: 


<y^^> = p = m/du- 

(mass density) 

(18a) 

^"^ 2 ^ ~ P “ m/fudu 

(momentum density) 

(18b) 

<y^> = q = ^ /fu^du 

Ckinetic energy density) 

(18c) 




One can linearize the function f in the following form: 


f(u,r,t) = f^(u) [1 + h(u,r,t)J (19) 

The functional relation f^ is the equilibrium distribution given by 
the Maxwell distribution as 

f (u) = n (m/27Tk„T ) (-mu^/2k T ) (20) 

o o Bo 'Bo 

Therefore, the function h(u,r, t) is the perturbed quantity. If this 

expression for f is substituted into the governing Boltzman equation. 
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one obtains 


(|_ + u*V) h = (|^) (21) 

COLL 

Hence, the local density fluctuation is given by 

6p(v,t) *= itt/fQ (u)h (u,r , t) du (22) 

Therefore, to determine the density fluctuation one needs to 
determine the expression for the perturbation function h. This 
procedure is a long, tedious mathematical journey. Since it is 
not the purpose of this work to rederive the existing theory, it 
will not be given here. References 8,9 and the papers listed in 
the Bibliography, can be examined for a more thorough explanation. 

Ill - EXPERIMENTAL CRITERIA 

In the discussion that follows, the specxf ications required of 
the experimental equipment are given. These requirements are actu- 
ally those needed to measure a Rayleigh-Brillouin spectrum, but odd- 
ly enough, most of these requirements (mainly the stability) are al- 
so needed to monitor density fluctuations which is the purpose of 
this investigation. In addition to these basic specifications, one 
major requirement of the frequency filtering device is also needed. 
This requirement will be discussed after the basic specifications are 
given. 

As in^a^-' light scattering techniques, the need for a reliable 
and stable monochromatic light source is obvious. Due to the nature 
of the Rayleigh-Brillouin phenomenon, the above requirements should 
be redefined and new specifications examined. These requirements can 
be obtained by considering the difference of Brillouin scattering 
over the various other techniques. 

In Rayleigh-Brillouin scattering, the entire spectrum is con- 
fined to a very small frequency range. One can visualize this if a 
simple calculation is performed for the spectrum of a typical sub- 
stance. If nitrogen is taken as an example, at T - 300 °K the shift 




becomes 97.07MHz and the halfwidth 34.2MHz. 

The next important factor is the low light level of the scattered 

radiation. This may seem to be a contradiction to what has been said 

earlier, but it is not. The reason being, in most light scattering 

density measurement methods a high power pulsed laser is used. (For 

example in Raman scattering techniques) . Here a relatively low power 

continuous wave laser is used. The CW laser was selected since one 

is interested in a continuous spectrum. Also its large scattering 

cross-section should compensate for the reduction in power. Through 

the course of this experiment (and through some literature re- 
8 9 

search) , ' it is discovered that either long running times, 15-30 
minutes, or repetitive scans were necessary to detect a spectrum. 
Therefore the larger Brillouin scattering cross-section did not com- 
pletely compensate the reduction in power. 

Keeping the above aspects of the Rayleigh-Brillouin scattering 
process in mind, some basic equipment requirements can be established. 

In regard to reliability, one cannot expect trouble free use of 
a highly complex instrument such as a laser. But after a period of 
debugging, it should operate with a minimum of problems providing 
the proper maintenance is performed. 

As for stability, there is actually two main types. These are 
energy output and frequency stability. In CW light scattering tech- 
niques, it is usually difficult to take into account variations in 
incident light energy. This involves some sort Of two channel cor- 
relation. Therefore it is usgfiilly required that the energy doesn't 
fluctuate more than a few percent over the period of the experi- 
ment. In some applications where a fluctuating parameter is 
being monitored, the laser fluctuation can outweigh the physical 
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variation. 


Before the frequency stability requirement is discussed the 
area of laser line width should be explored. As can be seen by 
the calculations presented earlier, the laser line width should 
be less than 68.4MHz (twice the halfwidth). Typically this is 
not found in most lasers without the use of an etalon. Care must 
be taken when incorporating an etalon. This adds an additional 
area where the stability of the system can be affected. Returning 
to the question of laser line width, one could ask how narrow 
should the width be? This is difficult to estimate. A safe value 
of 30-40MIIZ is advised. 

Since long running times or repetitive scans are needed, a 
frequency drift would cause an artificial broadening of the bands. 
Therefore for accurate results, the drift should be kept at a mini- 
mum. For this particular experiment a stability of 50MHz/hr is 
within the limits of experimental error. 

The above stability requirements of the laser system are neces- 
sary, but stability must be maintained in the other components of the 
experiment as well. In high resolution techniques an interferometer 
is normally used. 

Since the operation of an interferometer requires an exact spacing 
of its mirror (on the order of X/2 or less) , it is very susceptible to 
thermal and mechanical oscillations. Therefore, the above instabili- 
ties must be controlled such that they produce drifts in frequency 
that are compatible with the laser system. 

Due to low light level, a means of enhancing the detection 
of the spectrum is necessary. There are two primary methods in sig- 
nal acquisition that can be used here. First long running times can 
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increase the quantity of light collected. A photon counting tube 
and system could be added to further aid the light collection. 

Secondly, a multi-scan analyzer with the ability to average out 
random signals could be utilized. The additional requirement 
needed to measure density fluctuations lies within the interferometer. 

To measure these fluctuations, the interferometer must act like a 
very narrow bandpass filter and select a single frequency within 
the Rayleigh-Brillouin spectrum. In order to accomplish this, the 
interferometer should be prebiased at a particular frequency by ap- 
plying a single voltage to the prezoelectric drive. 

IV - EXPERIMENTAL EQUIPMENT 

The experimental setup is seen schematically in Fig. 5 and pic- 
torially in Fig. 6. During this experiment, two laser systems v;ere 
used. The first system consisted of a Control Laser Model 554A argon- 
ion laser pumping a Coherent Model 590 dye laser. The unique property 
of most dye lasers is that its wavelength is variable. The wavelength 
selection is accomplished by changing the angular orientation of a 
birefringent filter. This filter is internally mounted (see Fig. 7) 
and can be moved by turning a micrometer type adjusting knob. Using 
all lines of the argon-ion laser, which amounted to about 8 watts, 
on the dye laser produced approximately 1.4 watts. However, the incorpor- 
ation of an intracavity etalon system (see Fig. 7) into the dye laser 
which is necessary to reduce the line width to the specified lOMHz, 
reduced the power to about 300-500mW. As it turned out, this inter- 
ferometer only reduced the power available. It did not reduce the 
line width to the specified bandwidth for any appreciable period of 
time. The stability was such that it was found useless for the pur- 
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pose at hand. 

The second laser system was just the argon-ion laser itself. 

The specifications called for approximately 3.5 watts single moded 
radiation at 514.5 namometars. This turned out to be an empty 
claim. To reduce its bandwidth, an intracavity etalon was installed. 
Since it is characteristic of argon-ion lasers to give off a large 
amount of heat, the laser utilizes a water jacket for cooling. 

Regardless of the source of illumination used, the laser light 
was directed toward the scattering volume by a Newport Research 
Corporation Model 670 beam steering mirror. The mirror is able to 
adjust the angular, horizontal and verticle position of the inci- 
dent light. The need for this precision mirror was justified by 
the extreme dependence of the experiment on an accurate alignment. 

To fulfill this requirement the beam had to be positioned exactly. 
This was sometimes difficult since it was often necessary to peak 
the laser to its optimum power level. By doing this, the spatial 
position shifted slightly. Therefore the mirror was needed to re- 
direct the beam to its rightful course. 

To investigate a pure substance a pressure cell v/as constructed. 
There was nothing unique about the cell. It was constructed of 
aluminum with a port on each face. To further reduce the effects of 
reflection, the cell was black anodized and the windows anti-reflec- 
tion coated. The latter was not performed until much later in the 
investigation. The reason being it was not known how much the re- 
flections affected the background photon count. 

The light scattered from the material in question was collected 
by a 5cm diameter, 5cm focal length lens. To produce a point image 
of the scattering volume, the lens was positioned two focal lengths 
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from the scattering volume. At this point, an adjustable iris was 
placed. This was to act as a spatial filter and limit the collec- 
ted Ixght to that coming from the scattering volume. Another 5cm 
lens was ‘placed 5cm from the iris, to produce a collimated beam of 
light. 

To frequency filter the scattered light, A Burleigh Model RC- 
110 interferometer was used. As a means of frequency scanning, a 
piezoelectric drive was utilized. Since this model incorporated 
plane mirrors the free spectral range (FSR) could be varied. The 
FSR could be varied from IGHz to 15GHz by varing the mirror spacing 
from 15cm to 0.1cm. Because the parallelity of the mirrors is high- 
ly temperature dependent, the interferometer was kept in a thermally 
insulated, heated box. This box kept the interferometer at a con- 
stant temperature, 22+ ,05°C. 

As in all light scattering techniques a means of low light level 
detection is necessary. Here a RCA C31034-01 photon counting photo- 
multiplier tube was used. 

To reduce the background count and to electromagnetically shield 
the photomultiplier tube, the tube was housed in a Pacific Photo- 
metric thermoelectric ooler. Typically a noise of 300 counts per sec 
ond at a cathode voltage of 1800 volts was observed. 

As a means of signal processing several systems were used. The 
system primarily utilized was a Princeton Applied Research (PAR) pho- 
ton counting system. It consisted of a model 1121 amplifer-discrim- 
inator and a model 1112 photon counter-processor. This system is 
capable of monitoring two input channels, by the use of a gated samp- 
ler. The two signals can then be either added or subtracted. The re 
suit is either displayed digitally or interfaced to a computer system 



During a short time span two other systems were used. These 
were Princeton Applied Research Model 4202 signal averager and 
Model 162 boxcar integrator. The first can accumulate data from 
repetitive scans. Thus the result is the average of all scans re- 
ceived by the tinit. This averaging process can be several types; 
normal, summation and exponential. 

The boxcar integrator provided similar results by a different 
method; after a series of scans are accumulated, the complete data 
is integrated to produce a single sweep where the random signals 
are integrated out. 

V - GAS SAMPLE PREPARATION AND TECHNIQUE 

To minimize the effects of dust and the presence of other 
gases the test chamber was evacuated and filled several times. A 
Veeco vacuum pressure transducer recorded a pressure of 20 microns 
absolute during the evacuation of the chamber. As an added pre- 
caution the gas sample was passed through a Matheson gas filter to 
eliminate dust particles. 

The cell was usually pressurized between 50-70 psia. This 
pressure was measured by a Heise mechanical pressure gauge. 

Before and during the experimental test runs the laser's band- 
width was monitored by a Jodon Model 1500 spectrum analyzer. 

In monitoring the density fluctuations, a modification in the 
usual operation of the interferometer was needed. To monitor the 
fluctuations one frequency must be selected. This can be accom- 
plished by presetting the piezoelectric drive at a particular voltgage. 
At this point the interferometer mirror spacing is fixed and there- 
fore it will only transmit light at one frequency. The value of 
this particular voltage can be obtained by simultaneously measuring 
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the spectral response of the laser light source, Fig. 8A, and the 

voltage applied to the interferometer, Fig. 8B. From this one can 

select any frequency by picking a particular voltage. The voltage 
*• 

selection depends on the region concerned. In the hydrodynamic 
region the obvious choice is the voltage at which the maximum in- 
tensity occurs. In the kinetic region a voltage should be selected 
which is free from the effects of the central component. A fair 
estimate of this can be obtained by calculating the appropriate 
shift. 

VI - FINDINGS 

The primary concern of this investigation was to assess the 
potential of Brillouin scattering as a diagnostic tool. As a first 
step toward this goal, a spectrum from a sample substance was sought. 
This would lay the ground work for later work by indicating the 
troublesome areas. 

Although this investigation could not produce a Brillouin spec- 
trum, the shortcomings of this method as a means of measuring in- 
stantaneous density fluctuations still surfaced. By analyzing the 
reasons for failure to detect the scattering process these short- 
comings become visible. Therefore a brief discussion of the problems 
encountered will be given first. Then the shortcomings will be dis- 
cussed in the conclusion. 

4 

As stated before, an earlier study attempted to measure density 
fluctuations in liquids. The results obtained from this investigation 
are given in Figs. 11 and 12. From these results the impact of the 
impurities present in the liquid could be seen. They produce a level 
of noise that makes the detection of the spectrum very difficult. 

Also, due to the inferior quality of the interferometer, this noise 
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could not be filtered out. It was for these reasons that new 
equipment was acquired. With this equipment a new study was under- 
taken. The results of this study are given below. 

In the first experimental setup the dye laser system (argon-ion 
and dye laser) , Burleigh interferometer and PAR photon counting 
system was utilized. These instruments were mounted on two work 
benches fitted with 1-inch thick aluminum table tops. It took seven- 
eight months to get the two lasers operating simultaneously. The 
problem here was the argon-ion power instability; and improper dve 
laser mirror and jet alignment. To correct these defects, a new 
argon-ion laser was obtained and the dye laser was shipped back for a 
proper alignment. 

As an initial check-out procedure of the dye laser, a prelimin- 
ary examination of the laser line width was undertaken. The line 
width was measured using the Burleigh interferometer. An example of 
this spectral display is seen in Fig. 9. Here the FSR is 3GHz. There- 
fore the laser line width is measured to be approximately 250 MHz. 

This IS not indicative of the best . value achievable. A bandwidth 
slightly under lOOMHz has been obtained. Even this value is far above 
the specifications set by the manufacturer. The reason for this will 
be discussed shortly. Even this "narrow" line width could not be main- 
tained for an appreciable time period. The line width had a tendency 
to broaden and drift in position. This drift would cause the number 
of orders to vary. Therefore it would vary the number of Rayleigh lines 
per Sweep. Also the laser line broadening would make the central peak 
appear wider than it should. Not only did the line width degenerate, 
but the intensity varied. Ascertaining these variations as a function 
of time proved to be difficult if not an impossible proposition. The 
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power meter used for this purpose had too large a response time to 

indicate the actual deviation. As far as the apparent frequency 

drift and spectral line broadening are concerned, it is not known 
** 

to what extent the temperature instability of the interferometer 
was responsible. Making the assumption that the insulated ther- 
mal box of the interferometer provided a better degree of stabil- 
ity, one must conclude that the time stability of the laser was 
at most ten minutes. 

The above instabilities were primarily due to thermal and mech- 
anical induced vibrations. The thermal effects present in the dye 
laser would allow the etalons to creep and thus destroy their 
parallelism. These temperature effects were the result of room temp- 
erature gradients. The causes of the mechanical vibrations were more 
difficult to detect. These were discovered to originate from the 
argon-ion's high voltage power transformer, pulsating cooling water, 
and the dye laser's dye pump. The most damaging of these was the 
♦ water supply and dye pump. The reason being the vibrations were 
I transmitted thorugh their respective supply lines which are connected 
directly to the lasers. A flexible tygon hose was connected in series 
with the dye supply line. This relieved some of the vibration but 
not all. What was really needed was a means of completely isolating 
the system from all these effects. This could be accomplished if the 
system was set on a vibration absorbing lab table. Due to their ex- 
orbitantly high cost one could not be obtained at that time. Hence, 
the experiments had to be run during the brief periods of "stability". 
These periods were either infrequent or not long enough to produce 
satisfactory results. 

The results that were obtained indicated a high background read- 
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ing. It was assumed either a misalignmert in- the optical system or a 
spatial laser drift caused the unusually hiah background count. Actixally 
these effects are strongly coupled and both probably contri- 
buted sumultaneously as well as separately. - laser drift could 
not be helped. It was a fact that one had to contend with. The 
high background count was thus the result of reflections and Mie 
scattering from impurities. Their effects could not be filtered out 
so readily by the interferometer since they are of the same wave- 
length and are orders of magnitude large in intensity. Further 
light baffels and traps were needed and incorporated. 

It was at this time a new experimental area was obtained. In 
this area was a device that appeared that would solve the majority 
of the problems encountered. This device was a 17,000 pound steel 
table that measured "'2 feet by 4 feet. It had the ability to isolate 
disturbances transmitted via the ground [since it could be supported 
by 16 inflated aircraft inner tubes] . 

The results that followed didn't appear any different than 
those obtained earlier. These results were still prone to laser in- 
stability and high background count. It was discovered the vibra- 
tions still prevailed. These vibrations were not damped out by the 
table because they were transmitted through the supply lines and not 
the ground. Also the temperature effects were still present. There 
was no means of eliminating them other than temperature controlling 
the system or the room. 

Since it appeared that the dye laser system could not be used 
because of its instability, the decision was made to use solely the 
argon-ion laser with its inter-cavity etalon. Initially the laser 
was not chosen for this purpose. Therefore no real consideration was 
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given to the etalon's performance. To obtain some information 
on its performance, an inquiry v;as made of the manufacturer. 

Their reply was not too favorable. They had no specifications 
as to its exact performance. Also they reluctantly divulged 
that due to its poor performance (stability and selectivity) they 
are considering removing it from sales. 

Since no good information was available, some measurements 
were made. The laser line width was measured to be between 125- 
150 MHz. An example of this is seen in Fig. 10. Here the band- 
width is measured to be 125MHz. This system was a little more 
stable than the dye laser system. Again it was too difficult to 
determine the degree of stability and only an estimated time of 
10-15 minutes could be assumed. Frequency stability was greatly 
affected by the large amounts of heat radiating from the laser. 
Thermal equilibrium could never really be reached due to the erratic 
flow rate of the cooling water. An attempt was made to eliminate 
the pulsations by installing a large water reservoir in series with 
the water supply and the laser head. Before the benefits of this 
system could be assessed, an accident occurred. Somehow the flow 
stopped to the head and the discharge tube overheated. Accompanied 
by this overheating was an increase in gas pressure and, therefore, 
a loss in power. Since there is no means of reducing this pressure, 
the laser had to be shipped back for repairs. This was only one of 
many journeys back to the manufacturer. In all there were at least 
six. Four of these resulted in the return of new lasers because the 
old ones couldn't be repaired. These trips usually lasted for three- 
four weeks. When the laser was in operation it had a tendency to 
lose power. This problem resulted from a thin film deposited in the 
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Brewster windows by the outgasing of the plastic, which protected 
these windows. With the use of a dry nitrogen purge system, the 
formation of this film was eliminated. 

The subsequent experimental tests were still plagued with the 
high background. Many methods to reduce this background were used. 
Some of them included light baffels, small aperture iris' and black- 
ended enclosures. The count couldn't be reduced significantly. 

Since the above procedure proved unsuccessful, other means of 
elimxnating this background was sought. Considering the background 
count as a random process, then a method of averaging would reveal 
only the signal present. This type of signal processing would how- 
ever defeat the purpose of this investigation. One has to remember 
that the purpose of this investigation was to measure the fluctu- 
ations of the density in a flow field. To this end, a boxcar in- 
tegrator and a multi-scan analyzer were acquired. These units were 
on loan and therefore the time they were available was limited. 

A combination of this limited time and the infrequency of the laser 
stability contributed to the failure of this attempt. 

VII - CONCLUSION 

Upon examination of the findings presented in the previous sec- 
tion, the shortcomings of Brillouin scattering as a diagnostic tool 
became visible. 

• First, due to the low light levels either long running times or 
repetitive scans are necessary. This is obviously unsatisfactory 
for measuring instantaneous density fluctuations. Next, the usual 
environment in fluid dynamical experiments is not void of contami- 
nants. These particles would produce an overwhelming response that 



would hide the actual signal. Also, the geometrical bodies and 
configurations usually used would make it impossible to eliminate 
light duQ to reflections. In addition to the above reasons, 
there is the expense and practicality. It would cost a consider- 
able amount to isolate the system vibrationally . This cost does 
not appear to be of practical use when other methods are available. 
Finally, there is a very important point that is not readily visi- 
ble from the results presented. Since the interferometer is un- 
stable in the scanning mode, there is no reason to expect it would 
be stable in the preset mode. The instability would cause a drift 
which in turn would allow the intensity to vary. This variation 
could be interpreted as a fluctuation in density which is erroneous 
Therefore with the equipment presently available, the use of Bril- 
louin scattering as a diagnostic tool is in doubt. 

In conclusion, it seems that Brillouin scattering should be 
left to those who wish to measure sound speeds , thermodynamic 
properties (ratio of specific heats, etc.), and transport proper- 
ties Cshear viscosity, bulk viscosity, and thermal conductivity) 
of fluids. 
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FIG. 4 BRILLOUIN SPECTRUM AS A 
FUNCTION OF y (FROM REF 4) 
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FIG. 5 SCHEMATIC VIEW OF THE EXPERIMENTAL 
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FIG. 9 SPECTRAL RESPONSE OF DYE 

LASER; FREE SPECTRAL RANGE 
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APPARENT SPECTRUM OBTAINED FROM WATER 
(TAKEN FROM REE 4) 




FIG.I2 SPECTRAL RESPONSE FROM ABSOLUTE ALCOHOL 
(TAKEN FROM REE 4) 
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